Massive fields tend to form highly oscillating 
self-similarly expanding shells 



Gyula FodoiQ and Istvan Racjj] 

MTA KFKI, Reszecske- es Magfizikai Kutatdintezet 
H-1121 Budapest, Konkoly Thege Miklos lit 29-33. 
Hungary 

(Dated: February 1, 2008) 

The time evolution of self-interacting spherically symmetric scalar fields in Minkowski spacetime is investi- 
gated based on the use of Green's theorem. It is shown that a massive Klein-Gordon field can be characterized 
by the formation of certain expanding shell structures where all the shells are built up by very high frequency 
oscillations. This oscillation is found to be modulated by the product of a simple time decaying factor of the 
form t -3 / 2 and of an essentially self-similar expansion. Apart from this self-similar expansion the developed 
shell structure is preserved by the evolution. In particular, the energy transported by each shell appears to 
be time independent. 
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INTRODUCTION 



Because of its importance in black hole physics, the late 
time evolution of various fields on a fixed spherically sym- 
metric static asymptotically flat background spacetime is 
of obvious interest. The physical mechanism by which a 
scalar field in such an asymptotically flat spacetime is ra- 
diated away has been extensively investigated by means 
of analytical and numerical techniques. First, the case 
of massless fields was considered PJ 0, 0, 0, 0| later the 
evolution of self-interacting (massive) fields was also in- 
vestigated [HQ] IMS' I n both cases it appeared that the 
late time evolution of scalar fields is dominated by an in- 
verse power-law behavior. In particular, it has been found 
(see, e.g., Q) that at a fixed radius, in the intermediate 
asymptotic region, each multipole moment of a self- 
interacting scalar field $ with a mass parameter m evolves 
according to the oscillatory inverse power-law behavior 
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It is a common feature of all of the previous investiga- 
tions that the late time behavior of the scalar field is mon- 
itored only through the investigation of the field variable 
at a fixed radius. This might explain why the self-similar 
part of the decay rate, described in detail below, was not 
noticed in either of these considerations. 

There has been only very limited attention paid to the 
study of the behavior of massive fields at null infinity. 
Among the very few relevant investigations the most im- 
portant one concerning this problem is due to Winicour. 
He showed that the time evolution of massive fields in 
case of initial data of compact support necessarily yields 
0(l/r°°) asymptotic behavior at null infinity ^jj. 
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Having all the above results, it might be somewhat sur- 
prising that this paper is dedicated to the investigation of 
the dynamical properties of a massive spherically symmet- 
ric Klein-Gordon (KG) field on the Minkowski spacetime. 
In fact the necessity of studying this problem arose during 
the examination of the evolution of a more complicated 
dynamical system. We started by investigating the evo- 
lution of "excited magnetic monopoles". More precisely, 
in [Tl| the time evolution of spherically symmetric Yang- 
Mills-Higgs systems is numerically considered on a fixed 
Minkowski background spacetime. To have a computa- 
tional grid covering the full physical spacetime, the tech- 
niques of conformal compactification, along with the hy- 
perboloidal initial value problem, were used. In this way it 
was possible to study the asymptotic behavior of the fields 
close to future null infinity for considerably long physical 
time intervals. The numerical simulations indicated that 
when massive fields are included certain expanding shell 
structures form, where all the shells are built up by very 
high frequency oscillations. Since the maximal frequency 
of these oscillations is growing in (both the physical and 
nonphysical) time, the associated wavelengths reach the 
size of any fixed equal distance grid spacing specified in 
advance in the calculations. This may be expected to yield 
a significant contamination of the numerical data. There- 
fore, after only the first dozens of runs, we suspected that 
the formation of the above mentioned high frequency oscil- 
lations and expanding shell structures was due to numeri- 
cal noise. Nevertheless, convergence tests, along with the 
monitoring of the energy conservation, made it apparent 
that we have to take this as a physical phenomenon. The 
next natural question was whether the existence of these 
structures is due to the nonlinearities associated with a 
coupled Yang-Mills-Higgs system or whether the same ef- 
fect may also occur in the case of a simple linear field as 
well. This issue is of obvious physical importance, because 
far from the center of symmetry the field values are ex- 
pected to differ only slightly from that of the static "back- 
ground" monopole solution; moreover, the linearized equa- 
tions, relevant for the asymptotic region, take the form of 
a pair of uncoupled massive and massless KG fields. Our 
aim in this paper, in addition to presenting the main fea- 
tures of this phenomenon, is to demonstrate that they are 
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not the nonlinear effects that are responsible for the ap- 
pearance of expanding highly oscillating shell structures, 
since the same features already appear in the evolution of 
massive scalar fields. 

To simplify the physical setting, in this paper only the 
evolution of a spherical symmetric massive KG field is 
investigated. It is shown that for essentially arbitrary 
initial data with compact support the evolution can be 
characterized by the formation of expanding shells built 
up by very high frequency oscillations. As the time passes 
the maximal physical frequency of the oscillations forming 
the outer shells increases without any upper bound and 
thereby more and more shells become visible. There is 
also an argument presented explaining the qualitative and 
quantitative properties of the underlying physical process. 
We found that the exact time evolution of an initial data 
function with compact support can always be approxi- 
mated by a simple expression to a very high precision in 
a considerably large portion of the Cauchy development. 
It turns out that the evolution of such a massive KG field 
yields high frequency oscillations modulated by the prod- 
uct of a simple time decaying factor of the form t~ 3 ^ 2 and 
of an essentially self-similar expansion. Despite the mas- 
sive character of the KG field, the "edge" of this self-similar 
expansion moves with the velocity of light. This behavior 
of the KG field has also been justified for the full (non- 
approximate) description by numerical integration based 
on Green's theorem. As far as the authors know, the ap- 
pearance of this scaled down self-similar feature of the 
modulation of the high frequency mode has not been no- 
ticed yet, and it seems to be of physical interest in its 
own right. The developed shell structures are found to be 
stable, i.e., they seem to have their own lasting individual 
existence. In particular, the total energy which can be 
associated with each of them and which is transported by 
them is apparently conserved during the evolution despite 
the underlying expansion. 

The plan of this paper is as follows. In Sec. El we briefly 
describe the investigated physical system. In Sec. IHII the 
time evolution of the deformation of a vacuum configu- 
ration is considered. The exact and approximated field 
values are determined by making use of Green's theorem 
and a suitable approximation proc ess. The correspond- 
ing analysis is outlined in Sec. IIVI for temporarily static 
initial data. In Sec, fvl further characterization of the ob- 
served expanding shell structures is presented in terms of 
the energy density and energy current density profiles. We 
conclude in Sec. ED with a brief summary of our results 
along with some of their implications. 



II. PRELIMINARIES 



This section is to recall some of the basic notions and 
results we shall use in considering the evolution of a mas- 
sive KG field $ satisfying 



□ $ + m^$ = 



(HI) 



on a fixed Minkowski background. Given a point source 
at a point with Minkowski coordinates x' a , the retarded 
Green's function yields the generated field value at the 



point x as 

G(x\x) — —H(x° — x'°) 

ZTY 



S{\) - -^=H(X)Ji(mVX) 



2a/A 

(II.2) 

where H is the Heaviside function, 8 is the Dirac delta 
function, J\ is the Bessel function of the first kind of order 
one, and A is the square of the Lorentz distance separating 
the two points, i.e. 

A = (x a — x' a )(x a — x'a). (II.3) 

Using the Bessel function relationships Ji(£) = — Jj"A)(£) 
and Jo(0) = 1, the Green's function can also be written 
in the alternative form 

G (x;x') = -^H(x° - x' )-^ [jf(A)J„(mVA)] . (IL4) 

Given the field $ and its normal derivative n f d f $ on a 
spacelike hypersurface E, by virtue of Green's theorem, in 
the causal future of E the field value can be expressed as 

<£>(*) = jdE[G(x;x')n^<f(x')~<f(x')n^G(x;x')}, 

(11.5) 

where n e denotes the future pointed unit normal vector 
field on E. Whenever E is chosen to be a t = const, flat 
hypersurface of Minkowski spacetime the operator ri gJW 
is simply the partial derivative d t i . 

The first term in Eq. l|II.5|l gives the contribution 
yielded by the excitation (dt^)s of a vacuum initial data 
function with <E>|s = 0, while the second term can repre- 
sent the evolution of a temporarily static configuration, 
i.e., whenever (d t $)i: = but $|e is nonvanishing. As we 
will see in the following sections it is easier to consider the 
evolution of these two special types of initial data settings 
separately. Note that because of the linearity of the sys- 
tem the evolution of a general initial data specification is 
simply yielded by the superposition of these two types. 

Whenever both the background and the KG field are 
spherically symmetric it is advantageous to use spherical 
coordinates (t,r,i},ip). Then the line element reads 



dr - 

while Eq. Ijll.lfl takes the form 
d 2 $ d 2 § 2 <9$ 



d-ff 2 + sin 2 6 dip 2 ) 



dt 2 



dr 2 



dr 



+ m $ = 



(II.6) 



(II.7) 



for the field variable $ = «3?(t, r) 
density of the field is 



The associated energy 



<9£ 

dt 



+ 



9$ 
dr 



while the outgoing energy current density is 

'd<f>\ /<9$ N 
~ffi) \~dr 



(II. 



(II.9) 



In analyzing the behavior of the spherical field $, one 
can assume, without loss of generality, that the observer 
lies on the axis of rotation associated with the spheri- 
cal coordinate system, with coordinates i° = (t, r, 0, 0). 
Denoting the coordinates of t he source point as x' a — 
(t! ,r' ,ip'), the relation dll.'it takes the form 



A = (t — t') 2 — r' 2 — r 2 + 2r'r cos 



(11.10) 
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III. DEFORMATION OF A VACUUM 
CONFIGURATION 



We start off by choosing the t = hypersurface as our 
initial data surface Eo and moreover assume that $|e — 
but (d t $)\s = $o(r) where $ : [0, oo) — » R is a suffi- 
ciently regular function of r. To get the field values we 
have to evaluate only the first term in Eq. 111.51 . By mak- 
ing use of Eq. III.4I at a point of the axis of rotation with 
coordinates x a = (t, r, 0, 0) with t > 0, we get 

OO 7T 27T 

$(t,r) = j dr j M 1 j dtp'$ a (r')G(t,r;r' ',<p' ,#')r' 2 sm$' 



OO 7T 

r dr'r' 2 $ (r') / di?'-^ [i?(A) Jo(mVA)] sintf'. 
o 

("I 1) 

By virtue of Eq. I^TTTTTTfl . ^ = 

- 2r ; ram^ af 7 holds. This makes it 
the v integral, which provides 

OO 

*(t,r) = Yr\ dr ' r '*°( r ') [ff(Ao)Jo(mV^) 
o 

-H(A T )J (mVA^)| , (III.2) 



where Ao and \ n are the values of A at #' = and i?' = 7T, 
i.e., 

(III.3) 



A = t 2 -(r-r') 2 



(r + i 



(III.4) 



In general, the field value given by the expression (jlll . 2tl 
can be evaluated only numerically. This evaluation, how- 
ever, can be done to a high precision very efficiently, for 
example, by using the numerical integration package of 
the GNU Scientific Library Q. 

In the concrete examples described in detail in this pa- 
per the initial data functions are chosen to possess the 
form 



cexp 



Trrp-j , if r 6 [a- b, a + b]; 
0, otherwise, 



(III.5) 

for r > with center at r = a(> 0) and with width 2b > 0, 
which is a smooth function with compact support. In par- 
ticular, in this section the graphs presented refer to the 
evolution of an initial data function $ = f(r) with a = 4, 
b = 3, c = 1, and d — 10. This initial data function cor- 
responds to "hitting" the vacuum configuration between 
two concentric shells at r = 1 and r = 7 with a bell shape 
distribution. (The energy density and the energy current 
density profiles associated with the same initial configu- 
ration are shown by Figs. Q 03 and El below.) The field 
value 3> associated with m = 1 at the time level surfaces 
t = 100, t = 1000, and t = 10000 are as shown in Fig. 
U The final segment of the last two time slices is plotted 
in Fig. El In order to facilitate comparison of the differ- 
ent time slices, to compensate the fast expansion the field 



value is plotted as a function of the self-similarity vari- 
able p = j. The striking similarity of the structures at 
different time slices also holds for other choices of initial 
data functions, although the exact shape and appearance 
of the shell contours may differ significantly. 




FIG. 1: The time evolution of a KG field with initial data 
*|so = an< ^ (^**)|e = is shown. In order to compensate 
the fast expansion and the decay in time, instead of the field 
amplitude <& the combination t 3 / 2 $ is plotted on the three 
successive time slices t = 100, t = 1000, and t = 10000 against 
the variable p = j- Notice that the value p = 1 labels the point 
sitting at the light cone emanating from the center at t = 0. 
The graphs - for the definition of see 1111.171 and 

181 - are also shown by the dashed curves. On the early 
time slice t = 100 the final segment does not fit the contours 
; however, at later times where conditions lIII.ll and lIII.2l 
are satisfied in larger and larger portions of the time slices, the 
exact oscillations get closer and closer to the contours. 



Next, we describe a method based on an approximation 
which explains the appearance of these stable structures 
and provides a quantitative account of the underlying phe- 
nomenon to a high precision. It is surprising by itself that 
under the conditions introduced below the integrand of 
Eq. l|HI.2ll can be very precisely approximated by a func- 
tion so that the integral dill . 2tl takes a very simple form. 

We denote by f' the radius of the smallest sphere 
B(0, r') centered at the origin so that B(0, r) contains 
the compact support of the specified initial data on Eo. 

Condition III.l ; We shall say that a point p, with co- 
ordinates (t,r), of the time slice Et is in the region of the 
approximation if for a given e < 1 the relation 



< e(t-r) 



(III.6) 



holds. 
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FIG. 2: The final segments of the plots of Fig. Q] for p > 0.9 
on the time slices t = 1000 and t = 10000. It is remarkable to 
what extent the high frequency oscillations fit to the contours 
zt*^, even in the region far away from the origin on the time 
slice t = 10000. 



Notice that on a fixed time slice the smaller the value of 
e the more precise the approximated field value, although 
the domain of the approximation gets smaller. The nu- 
merical results indicate that even for the relatively large 
e = 1/10 there is surprisingly good agreement between the 
exact and approximated field values. For an arbitrarily 
small fixed value of e the domain where the approxima- 
tion breaks down is t — r /e < r < t + r' . The size of 
this region, independently of time, is f'(l + 1/e), while 
the domain of the approximation < t — f'/e grows lin- 
early with t. Thereby one may think of condition IIII.il 
as holding asymptotically almost everywhere inside the 
causal future of the compact support of the initial data. 

In the domain where cond it ion IHI.ll holds. both Ao and 
\ir are positive, so then H(\q) = //(A,r) = 1. In addition, 
we have that 



t 2 -r 2 (t-r)(t + r) (t-r) 
and similarly 



< e 



r r 
< t. 7— < e. 



f 2 - r 2 (t-r) t 



(III.7) 



(III.; 



Then, by making use of these relations, along with the 
notation 

\ (III.9) 



lo = Vt 2 - r 2 , 
Ao and can be approximated as 
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A?, 



lo 1 + 
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(111.10) 

(111.11) 



Moreover, for z 3> 1 the Bessel function Jo (z) can also be 
very closely approximated as 



(ni.i2) 



By virtue of Eqs. l|III.lfl|l and JIII.11L whenever condition 
IIII. II holds, both the arguments of Jo in Eq. JIII.2L i.e., 
mvAo and m-^/A7, are much greater than 1 provided that 
the following condition holds. 



Condition III.2 



is satisfied. 



The condition 
mlo » 1 



(in. 13) 



This condition certainly holds in the region of approx- 
imation of condition IIII. II if t is chosen to be sufficiently 
large, i.e., for t > i. 



For large values of z the function y ^ sin [z + j) has 
a high frequency oscillation with a negligible change in 
the amplitude. Therefore we shall further simplify the 
approximation JIII. 1 2t of Jo(-z) applied for the arguments 
z — m\/Ao and z = ra\f\^ as 



Jo («) 



■ sin z + 



(111.14) 



Hence, whenever conditions IIII. II and IIII. 21 are satisfied, 
the substitution of Eq. ljIII.14fl into Eq. HII.2(l yields 

r'$ (r') 



$(i,r) 



1 

27 



ivm lo 



(m\f\^ + ^jj dr. 

(111.15) 



Then by making use of Eqs. ijHI.lOll . 1III.1H and keep- 
ing the leading order term, a straightforward calculation 
results in 



$(t,r) 



2 / , 7T 

— cos ml + — 

■Km lo V 4 



r'<to(r') sin 



mr 
lo-' 



dr. (111.16) 



Finally, by introducing p = j, the self-similarity variable, 
Eq. l|HI.16|l can be recast into the form 



where 

*6» = 



r + T 

4 



(111.17) 



miv 
r $ (r) sin 



mpr 



dr' . 



(111.18) 



Notice that the integral term of Eq. flll.181 is, in fact, the 
real part of the Fourier transform of the function r$o(r). 
The cosine term of Eq. 1111.17^ gives a high frequency os- 
cillation that is completely independent of the specified 
initial data function, i.e., of $o(r). The amplitude of this 
high frequency oscillation is modulated by the rest of the 
expression. There is an overall factor t~ 3 ^ 2 scaling down 
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the modulation in time. However, the term (p) de- 
pends on t and r only in the combination p = r/t, so this 
term represents the self-similar outward expansion of the 
modulation. 

For simple polynomial choices of the initial data func- 
tion <£> , the explicit form of vE' c j o (p) can always be de- 
termined. In particular, for a centered step function, i.e., 
with $ taking the value c G R + for < r < r b and being 
zero elsewhere, 





(l-p 2 )t 


mn 


m 2 p z 



mrtp 



mr b p 



mr b p 

(111.19) 



It is remarkable to what extent the high frequency os- 
cillations fit the contours ±3' < j o (see Figs. [Jand[2jl. We 
would like to emphasize, however, that not only this over- 
all behavior can be described properly by the approxi- 
mation method outlined above. Figure 3 shows the exact 
field value <3>, along with its approximated value <fr a , where 
the values of <3> and $ a are determined by making use of 
the relations <III.2fl and {III.17I . respectively. Since there 
is an almost exact coincidence between the two relevant 
curves, only the final segment with p > 0.93 is shown in 
Fig.Elfor the time slice t = 1000. 




FIG. 3: The exact value field value $ along with its value 
<&a as it is approximated by the relation MI. 171 are shown for 
p > 0.93 and for the time slice t = 1000. 

We would like to emphasize that for a fixed radius r 
and for large enough values of t we have that (p) £s 

I oo 

**„(°) = V 2 ? J r' 2 $o(r')dr' and lo « t, so that for a 
v * o 

scalar monopole Eq. l|III.17|l restores the inverse power- 
law behavior of Eq. Ijl.lfl relevant for I — 0. 

One might have the impression that the apparent in- 
crease of the maximal frequency of the oscillations on the 
above figures is merely a coordinate effect, i.e., it is a sim- 
ple consequence of the use of the self-similarity variable 
p = r/t, and one might suspect that it will not occur if 



the plots are made against the radial coordinate r. Nev- 
ertheless, the following graphs (see Fig. show that the 
increase of the maximal frequency along the t = const. 
slices is a true physical effect. 
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FIG. 4: Equal length (Ar = 60) final segments of the true 
(nonapproximated) oscillations of i -3 / 2 <!> are shown at the time 
slices t = 1000, t = 10000, and t = 100000. The associated 
graphs demonstrate that the increase of the maximal frequency 
of the oscillations does really occur, even though the plots are 
made against the r coordinate representing real physical dis- 
tances. 

A simple explanation of this phenomenon can be given 
based on the use of the approximation l|III.17|l as follows. 
Since the term responsible for the high frequency oscilla- 
tions is cos(mlo(t,r) + tt/4) the frequencies ujt and uo r - 
relevant for the r — const, lines and for the t = const. 
hypersurfaces - can be read off from the approximation of 
the term mlo in the neighborhood of a point with coordi- 
nates (t, r) as 

ml (t + At, r + Ar) w ml (t, r) + u t At - co r Ar, (III. 20) 
where the frequencies are 

t 



Ut ~ mdtlo = m- 



= —mdflo — m 



(111.21) 
(111.22) 



respectively. It follows from these expressions that as r 
tends to the value of t both u)t and u r grow without any 
upper bound, although their ratio tends to 1. This, in par- 
ticular, implies that the self-similar part of the modulation 
expands (at least at the "outer edge") with the velocity of 
light. The frequency of the oscillation at a fixed distance 
Ar =t — r from the "outer edge" (t — r) grows as vt, as 
it can be seen from the approximation 

u>r ~ w t = 7 =J '- : w , 1 Vt (III. 23) 



1 Ar^/t + r 
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relevant for Ar <C t. This, in turn, implies - in accor- 
dance with Fig. that more and more oscillations will 
be associated with the final segments of length Ar as the 
time increases. 

It is also informative to consider what can be seen by 
a static observer, moving along an r = const, world-line 
with r 3> r . First, a very high frequency oscillation ar- 
rives at t = r — f with amplitude growing gradually from 
zero. In the domain where the approximation IIII.17I 
is valid, the amplitude of this oscillation becomes modu- 
lated by t~ 2 ^ (p), while p = r/t decreases from a value 
close to 1 to 0. By virtue of Eq. ljIII.21fl the frequency of 
the high frequency oscillation gradually decreases, settling 
down finally to the value uo t = m in accordance with Eq. 



IV. EVOLUTION OF A TEMPORARILY 
STATIC CONFIGURATION 



Start now with a static initial data function, i.e., as- 
sume that (9 t <E>)|s = and <fr|s = $o(r) where <3> : 
[0, oo) — *• R is a sufficiently regular function of r with 
compact support. Now to get the associated field values 
we have to evaluate only the second term of Eq. llll.5fl . 
In particular, at a point of the axis of rotation with co- 
ordinates x a = (t,r, 0,0), with t > 0, we obtain from Eq. 

(SjJ 



*(t,r) = 



dr' J d-d' J dip' 
•$ (r') 

oo 

= ~J dr'r'$ (r') j M' 





-^pG{t,r;t',r',tp',-&') 



r sin it 
t'=o 



d 

where the relation 



S(X) - -^=H(X)Ji(mVx) 
2v A 



(IV.l) 



d dX { d\\~ 1 d t-f d 



dt' dt' \dfi' d<&' r'r sin & dW 



(IV.2) 



along with Eq. III.2I . has a lso been used. By evaluating 
the #' integral of Eq. flv"Tt we get 



*(f,r) = 



OO 

/ dr' r'$ (r') 



— — j=- Ji(mV A ) )=^Ji(mVA^) 

^ V v Ao v Avr 



(IV.3) 



For determination of the parts of the above integral con- 
taining the Dirac delta terms we used the relations 



7' OO 

5{Xo)i>{r')dr' = J 5(X )^(r')dr' + J 5(X )^(r')dr' 

r 

<5(A )V'(r')(a r 'Ao)" 1 dAo 

2 

<5(Ao)V(r')(a r 'A )- 1 dAo 
ip(r') \ 1 / V(r' 



-H(r-t) 
2 v ' \ r-r 



2\r 



2t 



[ip(r + t) + H(r - t)ip(r - t)], 



r' — r+t 

(IV.4) 



and 



/ 



-H(t- 



r' + r 



H(t-r)^(t-r), 



(IV.5) 



where tp : [0, oo) — > R is considered to be a sufficiently 
regular but otherwise arbitrary function. Thereby, we get 
from Eq. flV3> 



Att 



dr 



+ — [(r + t)$ (r + t) + (r- t)$ (|r - t\)] . 

(IV.6) 

This integral has been evaluated numerically and the rel- 
evant plots on the time slices t = 500 and t = 5000 are 
shown in Figs. 00 and Here the initial da ta fun ction 
was chosen to possess the analytic form of l|III.5fl . i.e., 
$° = f( r ), with parameters a = 0, 6 = 5, c = 1 and 
d = 15 corresponding to a bell shaped initial configura- 
tion centered at the origin with radius r = 5. 

In the domain where conditions !!!!. ll and lHI. 2l are valid, 
in particular, whenever t is sufficiently large, the nonin- 
tegral terms of Eq. l|IV.fi|l are zero, since both r + t and 
\r — t\ are larger than the radius ¥ of the initial data func- 
tion. Furthermore, the integral term in Eq. JIV.fill can 
be approximated by exactly the same procedure as that 
used in case of Eq. ^111.21 . To see this, note first that 
for z > 1 the Bessel function Ji(z) can be very closely 
approximated as 



t / \ / 2 / 7T 

Ji(z) « — \ — cos 2+- 

V TTZ V 4 



(IV.7) 



Then, in particular, whenever conditions IIII.ll and IIII.2I 
are satisfied by making use of Eqs. IIII.10I and Ijlll.llfl . 
along with the self-similarity variable p = j, the integral 
term, and thereby &(t, r) itself, can be approximated as 

$(t,r) w sin (m\/t 2 - r 2 + |) , (IV.8) 
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where 



, „ 2m 2n-2 

**o(p)=-\/— P 4 



■ / r' Q {r')s\n 



mpr 



dr' . 



(IV.9) 



FIG. 5: The time evolution of a KG field with initial data 
function 4>|e = ®° an d (dt^)\s = is plotted at the two 
successive time slices t = 500 and t = 5000 against the variable 
p = j . The graph s ±>I>.£ - for the definition of see Eqs. 
EVM and flvT9l - are also indicated by the dashed curves. 
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FIG. 6: The final segment of the second plot of Fig. |5] re- 
stricted to p > 0.9 together with an even smaller portion rele- 
vant for a later time slice t = 50000 with p > 0.995 are shown. 



Notice that Eq. 1IV.8II has the same type of structure as 
Eq. ljIII.17fl . The sine term of Eq. ijlV.Sfl again represents a 
high frequency oscillation that is completely independent 
of the specified initial data. The amplitude of this high 
frequency oscillation is modulated by the overall factor 
i -3//2 , scaling down the modulation in time, along with 
. Again the term depends on t and r o nly via 
the combination p = r/t and as in the case of Eq. HII.17|l 
this term gives a self-similar outward expansion of the 
modulation. 

V. THE APPEARANCE OF THE SHELLS IN 
TERMS OF THE ENERGY DENSITY 



As is obvious from the relations ill. Hi and 111.91 . to be 
able to evaluate the energy density and the energy current 
density expressions in addition to Eqs. 1III.2II and l|IV.6fl 
we also need to determine the derivatives of the field vari- 
able $ with respect to t and r, respectively. These deriva- 
tives, by straightforward calculation, although a bit more 
involved than the previous ones, for a generic initial data 
<&|e = $o(r) and (dt$)\i; = $o(r) can be shown to take 
the form 



ftSftr) = I$(t,r)-^ / ,'<!> (r'l 



A() V A,r 



dr 



+ - 



2r 



-'$o(r') 



dr' 



+ 7T \(r + t)$o{r + t) + (r - t)$ D (|r - t\) + (r + i)$' (r + 1) - \r - t\<5>' (\r - t\)] 
2r 



-i [«„(r + 1) - #„(|r -t\)]- 2-5 [(r + t)*„(r +t) + (r- t)$ a (\r - t\)] 



(V.l) 
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and 



+ _ / r '$ (r') 



2r 



■Ji(mVAo) j= Ji(mVXn) 



H(\ )(r-r') 
Ao 



dr 



J2(mVA ) J2(mVA^-J 



A. 



— IT ^ + ^ $ °( r + *) ~ ( r ~ *)*°0 - t\)] + — [(r + t)<S> (r +i)-(r- t)$ Q (|r - t\) 
+ $o(r + t) + (r + t)$' (r +t) + $„(|r - i|) + \r- t|$' (|r - *|)] - -$(t,r), 



(V.2) 



where stands for the r derivative of the function <& a (r). 
Then by making use of Eqs. ijll.Sfl and ijll.Qfl . along with 
Eqs. fllLlt , flV^ . and EH , we can evaluate 

the energy density and energy current density expressions. 
For the choice of initial data function used in Sec. IHII to 
produce Figs. U and |2j the various plots of the energy 
density on the time level surfaces t = 100, t = 1000, and 
t = 10000 are shown in Figs. and |U 

The corresponding energy current density profile on the 
time level surfaces t = 100, t = 1000 and t = 10000 is 
shown in Fig. El 

Similarly, for the choice of initial data function used in 
Sec. lIVI to produce Figs. Eland® the various plots of the 
energy density on the time level surfaces t = 500, t = 5000, 
and t = 50000 are shown in Figs. 1101 andllll 
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FIG. 7: The energy density profile relevant for t he excitation 
of the initially vacuum configuration of Sec. IIHI is shown for 
the time slices t = 100, t = 1000 and t = 10000. 



FIG. 8: The final segments of the energy density profile of 
Fig.[jjis shown for p > 0.9 and for the time slices t = 1000 and 
t = 10000. 



In order to represent more clearly the ratio of the total 
energy carried by individual shells in the spherically sym- 
metric configuration investigated, instead of the energy 
density e and the energy current density S, we plotted 
the quantities iirr 2 e and Anr 2 S representing the energy 
and energy current in thin spherical shells. Since we are 
using p — r/t as a radial coordinate on our plots we also 
have to multiply these values by j- = t 



E 



4:Tir edr - 



Aur tedp, 



(V.3) 



in order to be able to estimate the energy as the area 
under the plotted curves. 

Notice that the formation of self-similarly expanding 
shells is manifested strikingly by Figs. 171- 1111 in spite of 
the fact that no approximation has been used anywhere 
in the associated calculations. 



VI. CONCLUDING REMARKS 



In this paper the time evolution of initially concentrated 
spherically symmetric massive KG fields has been investi- 
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FIG. 9: The energy current density profile relevant for the 
excitation of the initially vacuum configuration of Sec. IIHI is 
shown for the time slices t = 100, t = 1000, and t = 10000. 
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FIG. 10: The energy density profile relevant for the excitation 
of the initially static configuration of Sec. IIVI is shown for the 
time slices t = 500 and t = 5000. 



gated. By means of concrete examples and by a suitable 
approximation method the following characteristic prop- 
erties have been found. There is an overall high frequency 
oscillation of the field value which is modulated by two 
factors. First, there is a time decaying factor of the form 
t~ 3 ^ 2 consistent with the total energy conservation. Sec- 
ond, the field amplitude is also modulated by a self-similar 
contour, i.e., apart from the scaling down t~ 3 ^ 2 the oscilla- 
tions expand essentially self-similarly. The associated self- 
similar contours were found to possess knots, whereby the 
oscillations are separated into individual shells. Since the 
overall exponent of the scaling down factor is independent 



FIG. 11: Two different final segments of the energy density 
profile of Fig. HJJ are shown for the time slice t = 5000 with 
p > 0.9 and for the time slice t = 50000 with p > 0.995. 



of the radial coordinate, the fraction of the energy carried 
by an individual expanding shell is also preserved during 
the entire evolution. The exact form of the shells, which is 
approximated to higher and higher precision in larger and 
larger domains of the associated domain of dependence, 
is determined essentially by the Fourier transform of the 
initial data. 

It is of obvious interest to know whether the results 
presented in this paper have any relevance for initial data 
with noncompact support. To be able to answer this ques- 
tion it is necessary to check whether the conditions applied 
throughout our analysis, i.e., conditions IIII.ll and IIII.2I 
can be satisfied at least in a generalized sense. We have 
investigated this issue briefly and we found the follow- 
ing. Suppose that we have a non-compact support initial 
data function, which, however, is focused sufficiently on 
the central region, i.e., it decays rapidly. Then we can 
associate a finite characteristic size with the support of 
such an initial data function. We expect that condition 
IIII.ll can be replaced by the requirement that this char- 
acteristic size should be small compared to e(t — r). If 
we have an initia l data function of this type and, in ad- 
dition, condition MI. 21 holds, basically the same type of 
approximations and plots can be produced for the field 
values and for the energy densities as was possible in the 
case of initial data with compact support. It might be 
interesting that for initial data functions $|e = $o(r) 
and (9t$)|s = 0, where <fr (r) was chosen to possess the 
form $ (r) = exp(— r 2 /d) (with d G there was no 

knot on the the associated contours ±^$ , while for the 
choice $ (r) = d/(d + r 8 ) (with d g R + ) the very same 
type of figures, with shells and oscillations, were obtained 
as in the case of initial data with compact support. In- 
vestigation of more general initial data specifications with 
noncompact support, as well as a careful adaptation of the 
approximation procedure applied throughout this paper, 
deserves further attention. 

It would also be important to know to what extent the 
results presented in this paper may have relevance in more 
complicated physical situations, for instance, in the case 
of nonlinear self-interacting fields, or whenever instead of 
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the Minkowski spacetime the evolution is investigated on a 
general asymptotically flat background spacetime. Note, 
however, that in the asymptotic region the field equa- 
tion relevant for a self-interacting possibly nonlinear scalar 
field in many cases is expected to reduce to the equation of 
a massive KG field; moreover, in that region the geometry 
also tends to the form of the Minkowski metric. Thereby 
it seems to be plausible to assume that the formation of 
shell structures, as well as the time decay rate and the self- 
similar behavior of the modulation, will probably occur in 
any asymptotically flat spacetime (regardless of whether 
the geometry is dynamical or not) and for any massive 
field that can be approximated by a massive KG field in 



the asymptotic region. 
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